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ABSTRACT 

We present the discovery of a wide (67 AU) substellar companion to 
the nearby (21 pc) young solar-metaUicity Ml dwarf CD— 35 2722, a mem- 
ber of the ;^100 Myr AB Doradus association. Two epochs of astrometry 
from the NICI Planet-Finding Campaign confirm that CD— 35 2722 B is phys- 
ically associated with the primary star. Near-IR spectra indicate a spectral 
type of L4±l with a moderately low surface gravity, making it one of the 
coolest young companions found to date. The absorption lines and near- 
IR continuum shape of CD— 35 2722 B agree especially well the dusty field 
L4.5 dwarf 2MASS J22244381 -0158521, while the near-IR colors and abso- 
lute magnitudes match those of the 5-Myr old L4 planetary- mass compan- 
ion, IRXS J160929. 1-210524 b. Overall, CD-35 2722 B appears to be an 
intermediate-age benchmark for L-dwarfs, with a less peaked H-hand contin- 
uum than the youngest objects and near-IR absorption lines comparable to field 
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objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find 
Teff =1700-1900 K and log(5') =4.5±0.5. The spectra also show that the ra- 
dial velocities of components A and B agree to within ±10 km/s, further con- 
firming their physical association. Using the age and bolometric luminosity of 
CD— 35 2722 B, we derive a mass of 31±8 Mjup from the Lyon/Dusty evolution- 
ary models. Altogether, young late-M to mid-L type companions appear to be 
overluminous for their near-IR spectral type compared to field objects, in contrast 
to the underluminosity of young late-L and early-T dwarfs. 



1. Introduction 



Currently, more than 440 extrasolar planets have been discovered by radial velocity and 



Butler et al. 
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Marcv et al. 


2008; 


Mavor et al. 


2009) 



(e.g. 

surveys are sensitive to Jupiter-mass planets at sep arations of > 10 AU around nearby stars 
25 pc) when they are young (< 100 Myrs) (e.g. 



Biller et al. 


2007; 


Lafreniere et al. 


2007 



Nielsen fc Close! l2010l ). So far, th ese searches have mostly yielded brown d\yarf or planetary 



mass objects at large separati ons ( Chauvin et al.ll2005al : iLuhman et al.ll2006l ; iLafreniere et al. 



20081 : iThalmann et al.l |2009| ). Spectroscopic studies of the atmosph eres of the brighter 



companions are reve aling their temperatures, gravities and masses ( iLavigne et al.l 12009 
Bonnefoy et al.ll2010[ ). However, even more exciting discoveries of planetary-mass objects 



at 5-120 AU separati ons are beginning to come forth ( iMarois et al.ll2008t iKalas et al.l 12008 



>ep 

Lagrange et al.ll2010D . Since these companions have separations that allow us to isolate their 
photons from those of the primary, we c an begin to charac terize the atmospheres of the 
companions discovered by direct imaging (IBowler et al.ll2010l ). 



Since December 2008, we have been conducting the Gemini NICI Planet-Finding Cam- 
paign to find and characterize planets through direct imaging (Liu et al. 2010). The Cam- 
paign is a three-year program at the Gemini-South 8.1-m Telescope to directly image and 
characterize gas-giant planets around 300 nearby stars. The NIC! instrument was specifi- 
cally designed for this purpose flFtaclas et al.ll2003t IChun et al.ll2008[ ). It is equipped with an 
85-element curvature adaptive optics system, a Lyot coronagraph, and a dual-camera system 
capable of simultaneous spectral difference imaging (IRacine et al.lll999l. SDF) on and off the 
1.6 /xm methane fe ature seen in cool (<1400 K) substellar objects (e.g. iBurrows et al.lll997l : 
Baraffe et al.ll2003l ). The two cameras are equipped with Aladdin II InSb detectors and have 
18.43" X 18.43" field of view and a plate scale of 18 mas/pixel (milli-arcseconds/pixel). NICI 
can also be used in a fixed Cassegrain Rotator mode for Angular Differential Imaging (ADI; 
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Liul l2004t iMarois et al.ll2005[ ). where the telescope image rotator is turned off, so that the 
slcy rotates on the detector while the telescope and instrument PSF (point spread function) 
artifacts remain fixed (See § 2.1.1). 

In an earlier paper we presented the discovery of PZ Tel B, one of few young substellar 
companions direct ly imaged at orbit al separations similar to those of giant planets in our 
own solar system ( iBiller et al.ll2010bl ). Here we present the discovery of CD— 35 2722 B, one 
of the coolest young (~ 100 Myr) close-separatio n (< 100 AU) brown dwarf companions to a 
star directly- imaged to date. [Torres et al.l ( l2008l ) identify CD— 35 2722 A as a member of the 
AB Doradus association, based on its common space motion and consistent galactic location 
with the other association members. They derive a kinematic distance estimate of ~ 24 pc 
(See § 3.1). The eponymous star of the AB Doradus associatio n itself appears to be coeval 



with t he Pleiades, which is thought be youn ger than 150 Myrs (IBarrado y Navascues et al. 



20041 : iBurke et al.1 |2004J : lOrtega et al.l 120071). and other studies conclud e that AB Dor is 



approximately 50 Myr or older (jZuckerman et al.ll2004l : IClose et al.l 120071 ). Thus, we adopt 
an age of 100±50 Myr for CD— 35 2722 A based on its membership in the AB Dor group. 



2. Observations and Data Reduction 



2.1. Gemini South/NICI imaging 

2.1.1. Campaign Observing Modes and Data Reduction 

We use two observing sequences for each target in the NICI Campaign to optimize our 
sensitivity to both methane- bearing and non-methane bea ring companions. Both sequences 
employ the ADI technique (jLiul l2004l : iMarois et al.ll2005l ). Exposure times per frame are 
kept around 1 minute so that the noise is background-limited, and the total time spent on 
detector readouts is small. Observing hour angles are chosen such that the rotation rate of 
the sky is less than 1 degree per minute to avoid excessive smearing during an individual 
exposure. The total sky rotation over the observing sequence is also required to be greater 
than 15 degrees^ to li mit flux loss by self-subtraction of companions at small separations 
(e.g. 



Biller et al.ll2009l ) 



To search for close-in methane-bearing planets, we combine the Angular and Spectral 
Difference Imaging methods into a single unified sequence ("ASDI"). In this observing mode, 
a 50-50 beam splitter in NICI divides the incoming light between the off (A = 1.578 /im) 
and on-methane (A = 1.652 /im) moderate-bandwidth (AA/A=4%) filters which pass the 
light into the two imaging cameras, henceforth designated the "blue" and "red" channels 
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respectively. The two cameras are read out simultaneously for each exposure. The red 
channel image is subtracted from the blue channel image to remove the nearly identical 
quasi-static speckles in the two channels. The resulting difference image has drastically 
improved sensitivity to methane-bearing objects, which are expected to be much fainter in 
the red channel (methane-absorbed) than in the blue channel. In the ASDI mode, we take 
45 1-minute frames. 

At larger separations (> 1.5"), our sensitivity is limited by throughput, not by speckle 
noise. Thus we also obtain an ADI sequence without the 50-50 beam splitter, sending all the 
light to the blue channel, in order to achieve maximum sensitivity. In this mode, we take 20 
1-minute frames using the standard H-hand filter, which is roughly four times wider than 
the 4% methane on/off filters. 

The initial image processing steps in the NICI Campaign data reduction pipeline were 
applied to all raw images. After bad pixel removal and flat-fielding, distortion corrections 
were applied. To determine the distortion, we acquired images of a rectangular pinhole 
grid placed upstream of the NICI AO system. Deviations of the imaged pinhole (spot) 
centers from a rectangular grid were used to measure the optical distortion of the blue 
and red channels. First, we fit a model rectangular grid of points to the spot centroids, 
allowing for rotation, change in grid spacing (different in x and y directions) and translation. 
Next, we found the 2nd order transformation in the coordinates of the observed spot centers 
that match the coordinates of the fitted model grid. The RMS in the difference between 
the distortion-corrected spot centroids and the rectangular model grid centroids was 0.1 
pixels. Distortion-corrected images were then constructed by replacing the values at the 
original coordinates with interpolated values at the transformed coordinates using cubic 
interpolation. 

In NICI data, the translucent focal plane mask causes arc-like features to appear near 
the half transmission radius (r = 0.32") of the Gaussian-profiled mask, which are dissimilar 
between the blue and red channels. Thus, for the ASDI and ADI datasets, we applied a filter 
to the images to remove arc-like PSF features centered on the primary star. For every pixel 
in the image, we computed the interpolated value at 100 points evenly spaced along an arc 
centered on the primary star and extending 30 pixels in azimuth around the chosen pixel. 
The median of these 100 values was then subtracted from the pixel value. This filtering 
process improves the quality of channel differencing during the ASDI and ADI reductions. 

An additional filter was applied to remove the stellar halo. Each individual image was 
smoothed with a Gaussian of 8 pixel FWHM and subtracted from the original (unsharp 
masking). Next, all the (unsharp masked) images were smoothed with a Gaussian of 2 
pixel FWHM. With this 2-step smoothing process, all spatial features that are significantly 
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different in size from the FWHM ~ 3 pixel NICI PSF were removed from the images. 



2.1.2. ADI and ASDI Observations 



CD— 35 2722 was first observed on UT January 16, 2009 in the standard campaign 
observing modes (ASDI and ADI). The NICI images revealed a bright companion candidate 
with separation ^ 3". For astrometric follow-up, we observed the system again on UT 
January 10, 2010 in the ASDI mode. Since the companion has separation > 1" and does not 
have methane, we reduced the blue and red channel data sets separately using a simple ADI 
data reduction method and then added the reduced images together. For each channel, all the 
images were registered using the primary star, which was detected at very high signal to noise 
through the partially transmissive focal plane mask. The images were median combined to 
build a PSF image. This image consisted only of quasi-static speckles, since any real objects 
in the dataset would appear at different PAs in the individual images (since the rotator is 
off) and would be removed by the median process. Next, the PSF image was matched to 
and subtracted from each of the 45 individual frames, fitting for the optimal translation 
and intensity normalization to match the speckle patterns between 0.65" and 0.78" from the 
star. This fitting region lies well outside the mask edge at 0.51". The final reduced image 
was constructed by rotating the individual subtracted images to place North up, registering 
them using the previously computed positions of the primary star, and median combining 
the stack. The PAs recorded in the FITS header of each frame give the orientation at the 
beginning of the exposures, but we calculated and used the PAs at mid-exposure to align our 
images. Finally, the reduced blue and red channel images were added together. We show a 
cutout from the final reduced image for the 2010 data in Figure [TJ 



2.1.3. JHKs Rotator- On Photometric Observations 

Upon detection of the very bright candidate at a separation of ~ 3", we obtained images 
in the JHKs filters on UT March 9, 2009 to measure its near infrared colors. The broad band 



filters in NICI are in the MKO (Mauna Kea Observatory) system ( ISimons fc Tokunagall2002 



Tokunaga &: Vaccall2005l ). The observations were obtained in the rotator-on mode. Nine 15- 
second exposures were taken, imaging the star simultaneously in both channels, with the H 
filter in the blue channel and the Kg filter in the red channel. Nine 15-second exposures 
with the J filter in the blue channel were taken immediately following these observations. 
For the reduction of the JHKs rotator-on data, we fiat-fielded and registered the frames on 
the stellar centroid and median combined the stack of images. Cutouts from the reduced 
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images of CD— 35 2722 in the J, H and X^-bands are shown in Figure [2l To determine 
the JHK band flux ratios of the companion relative to the primary, we first used simple 
aperture photometry (3 pixel radius) to measure the relative fluxes of the companion and 
the primary star, with the latter seen under the focal plane mask. 

To determine the transmission of NICI's coronagraphic mask, we imaged a pair of fleld 
stars, 4.2" apart on UT October 31, 2010. The stars, 2MASS J06180157-1412573 {H=QAQ 
mag) and 2MASS J06180179-1412599 (if=11.16 mag) were chosen so that they could be 
imaged at the same time without saturating either. Ten images were obtained in the H and 
i^-bands, with the brighter star under the mask, followed by 10 images with both stars off 
the mask. In the J-band, 4 images were taken in each conflguration. All images were taken 
with the Cassegrain rotator on, and thus there was no smearing during individual exposures. 
We then computed the ratio of the aperture (3 pixel radius) fluxes of the two stars, in both 
the on and off-mask conflgurations. The off-mask flux ratio divided by the on-mask flux ratio 
yields the mask attenuation. We determined that the central attenuation of the 0.32" mask 
is 6.37±0.06, 5.94±0.02 and 5.70±0.03 magnitudes in the J, H and K5-bands, respectively. 
The uncertainties are the RMS of the central attenuation measurements for the individual 
images in each fllter. 

These attenuation factors were applied to the raw flux ratio to arrive at the flnal JHKs 
flux ratios. We computed the uncertainty in the flux ratios by including the dispersion in 
the flux ratios over the data set and the uncertainties in the focal plane mask transmission. 
To convert the K^-band photometry of the companion to the -ftT - band fllter (a very small 



effect), we applied the transformations from [Stephens fc LeggettI ( 12004| ). Table [T] lists the 
final NICI photometry of CD-35 2722 B. 



2.2. IRTF/SpeX Imaging 

We also obtained images of CD— 35 2722 with the infrared guide camera of the facility 
spectrograph, SpeX, of the 3-m NASA Infrared Telescope Facility (IRTF) on UT February 
3, 2010 and UT September 24, 2010 in good seeing conditions. We obtained images in the 
JHK bands, using 5 dither positions for each filter. The IRTF/SpeX images were fiat- 
fielded, sky-subtracted, registered and stacked. To measure the fiux ratio of the primary and 
the companion we first removed the background fiux near the companion that comes from 
the primary. This was achieved by applying the same filter used in § 12.1.11 to remove arc- 
like features in the NICI PSF, except that medians were taken over arcs of length 90 pixels, 
instead of 30. Then fiux ratios between the primary and the companion were computed using 
apertures of 6 pixel radii. To assess the error in our measurements and any fiux loss from the 
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filtering process, a hundred simulated companion images were created by scaling down the 
primary to the flux level of the secondary and shifting them to the same separation, placed 
with uniform aziumuthal spacing, but at >20° difference in PA (>10 pixels away) compared 
to the actual companion. Photometry was done on the known positions of the simulated 
companions after the filtering process was applied. The RMS in the measurements was taken 
as the photometric error. We also computed the estimated flux loss from the mean of the 
differences between the known and measured photometry of the 100 simulated companions. 
These loss fractions were ~ 0.03 magnitudes, and we corrected for this loss in each band. 
Even though the two epochs of IRTF measurements agree to better than 0.02 magnitudes, to 
be conservative we adopt the errors from the simulations, and list the mean of the photometry 
from the two epochs in Table [TJ 

The primary star could be variable and thus use of the non-contemporaneous 2MASS 
photometry to derive the magnitudes of CD— 35 2722 B from our flux ratio measurements 
may not be appropriate. Thus we observed a standard star in JHK during our 2nd epoch 
IRTF observations on UT September 24, 2010. We measured JHK MKO magnitudes of 
7.85±0.03, 7.31±0.05 and 7.07±0.01 mag for CD-35 2722 A, which agree with the 2MASS- 
derived measurements to within la (see Table [T]). Moreover, other available data do not 
strongly indicate variability. Optical monitoring of CD— 35 2 722 A shows variability with an 



ampl itude of AV = 0.11-0.12 mag and a period of 1.7 days (jPojmanskill2003l : iMessina et al. 



2010l ). This is likely due to rotating spots on the surface of this young star, in which case 



the photometric amplitude should be smaller in the near-IR. 



2.3. Gemini North/NIFS Spectra 



Spectroscopic follow-up observations were obtained with NIFS flMcGregor et al.ll2003[ ) 
on the Gemini- North telescope. The target was observed with the J, H and K settings of 
NIFS, which cover 1.15-1.33 fim, 1.49-1.80 //m and 1.99-2.40 ^m, respectively, at spectral 
resolutions of AA/A = 5000-6000. Table [3] gives a log of the observations. The primary star 
was used as the AO natural guide star, providing significant correction of the atmospheric 
turbulence. Due to the southern declination of the target (DEC ~ -36°), the observations 
were carried out at high airmasses of 1.8-2.0. The average FWHM obtained was 0.15". 

The data were taken in a 'ABBA' nodding sequence, with the object dithered by 1.5" 
between two positions in the field of view. The data were reduced with the Gemini NIFS 
IRAF pipeline, following the standard procedures to produce a wavelength and distortion- 
corrected data cube. The sky background was removed from A frames by subtracting the 
negative of the nearest B frame, and vice versa. Figure [3] shows CD— 35 2722 B in an A— B 
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NIFS image. 

Spectra of the object were then extracted from each cube with the nf extract routine, 
using an aperture of 0.5". The wavelength solution for the 1-d spectra was checked against 
night sky lines (by re-running the reduction procedure without the sky subtraction), and 
re-interpolated where necessary. We found small drifts in the wavelength solution between 
nights on the order of ^ 2 A. 

Telluric standard stars of spectral type AOV were observed at similar airmass to the 
science observations and were reduced in the s ame way. Stellar fe atures (hydrogen lines) 



were removed with the IDL package spextool fjCushing et al.ll2004l ). The resulting telluric 



spectrum was then used to correct each individual science spectrum before combining them. 

The individual science spectra were scaled and combined using a 3-sigma clipped mean, 
with the rms per wavelength bin used to estimate the noise in the resulting combined spec- 
trum. The combined spectrum was then multiplied by the filter curves for the Gemini JHK 
filters, while the total flux in each band was set to the JHK photometry measured by NICI. 
The calibrated JHK spectra are shown in Figure |H 



3. Results 

3.1. Astrometric Confirmation of CD-35 2722 B 

Two objects were detected in the CD— 35 2722 field at separations of 3.172ib0.005" and 
5.312zb0.011" from CD— 35 2722 A on UT January 16, 2009 and again at separations of 
3.137±0.005" and 5.316±0.011" on UT January 10, 2010. Since the coronagraphic mask is 
not completely opaque, we should be able to measure the position of both the primary star 
and the candidate to better than a tenth of a pixel, given the signal to noise of the detections 
and the accuracy of our centroiding algorithm. However, our astrometry is limited by the 
uncertainty in our measurement of the image distortion. Comparison of two epochs of dense 
star field images shows the astrometric discrepancy between epochs was ~ 0.3 pixel at 3" 
and ~ 0.6 pixel at 6" and we incorporate these errors in our astrometry in Table [21 

To check whether these astrometric measurements are consistent with the motions of 
background objects, we need the proper motion and the distance to the primary. The 
proper motio ns for CD— 35 2722 A in RA and DEC are — 5.6±0.9 and — 56.6±0.9 mas/year. 



respectively ( IZacharias et al.ll2009l ). 



Using the CAPScam/duPont telescope at Las Campanas Observatory (IBoss et al.ll2009l ) 



we obtained 3 epochs of parallax observations separated by 8 months (April 2010, November 



- 10 - 



2010, December 2010). This yielded a trigonometric parallax of 47±3 mas (21.3±1.4 pc). 
Three is the minimum number of observations required to measure a parallax. Thus it is 
difficult to properly assess the uncertainty in our measurement. So we estimated the error 
using a Monte Carlo approach, assuming a conservative 1.5 mas precision on each epoch and 
fitting the parallax and proper motion of 10,000 Monte Carlo realizations of the data. Further 
observations are in progress and a definitive value will be published in Anglada-Escude et 
al. 2010 (in preparation). 

The mean distance to the AB Doradus group nucleus h as been refined in a new investi- 
gation of the moving groups within 100 pc of the Earth (see iMama j eki |20 1 Ol ) . The kinematic 
analysis therein employs th e space motions of only the group members with Hipparcos par- 
allaxes fjvan LeeuwenI 120071 ) and the highest quality radial velocity measurements to find 
a distance of 20.6±2.4 pc for 8 "nuclear" members. Comparing the AB Dor group space 
motions and the CD— 35 2722 A proper motion, we see that the star is very likely also a 
"nuclear" AB Doradus member and thus at the same distance. We note that this distance 
estimate is in agreement with our parallax measurement. The predicted radial velocity from 
the best-fit space motions is 31.5 ±1.4 km/s. This is completely consistent with the published 
measurement of 31.4±0.4 km/s ( iTorres et al.l 120061 ). 



In Figure O we compare the motions of the two objects in the CD— 35 2722 field to the 
expected motion of a background object given the first epoch astrometry, the known proper 
motion and estimated distance to CD— 35 2722 A. The change in the PA of CD— 35 2722 B 
is significantly different (at the 3a level) from that expected of a background object. This 
indicates that CD— 35 2722 B is physically associated with the primary, although about 
1.6 AU/year of orbital motion is required to explain the relative motion with CD— 35 2722 A. 
In comparison, the second object in the field at ^ 5.3" separation moves as expected for a 
background object, thus confirming the joint accuracy of the proper motion and parallax of 
CD— 35 2722 A and our NICI astrometric calibration. (See also Biller et al. 2010.) 

Since the viewing angle and eccentricity of companion's orbit are unknown, w e estimate 
its semi- major axis by applying a statistical correction from iDupuy Sz Liul (120101 ). Assum- 
ing a uniform distribution of orbital eccentricities from to 1, we multiply the projected 
physical separation by l.l^Qgg to estimate the true semi-major axis. Applying this factor 
to CD— 35 2722 B's measured separation, we obtain an estimated orbital semi-major axis of 
83lq o AU. We derived a mass of O.4M0 for CD— 35 2722 A using its age and absolute H mag 



with 



Siess et al.l (120001 ) evolutionary models. Using this mass with Kepler's Third Law we 



obtain a period of 1195_7]^q years for component B. The observed motion of the companion 
(1.6 AU over 1 year) is plausible for a bound object. 
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3.2. Near-IR Colors and Metallicity 



To measure colors of CD— 35 2722 A, we convert our relative flux ratio measurements 
(§ 2.1.3 and § 2.2) into apparent magnitudes for CD— 35 2722 B using the 2MASS photom- 
etry. The JHK-ha.nd IRTF/SpeX data yielded more accurate flux ratios than the NICI 
observations, because the IRTF photometry did not incur the additional uncertainty of the 
NICI focal plane mask transmission; thus, we use the IRTF-derived results for the photom- 
etry of CD— 35 2722 B. However, the NICI images provide far better astrometry than the 
IRTF data. Table [1] summarizes the resulting magnitudes for CD— 35 2722 B. 

Figure E] compares the near-IR colors of CD— 35 2722 B with other stellar and substellar 
objects, including old (high-gravity) field ultracool dwarfs and young (lower-gravity) objects 
found in the field and as companions. The JHK colors of CD— 35 2722 B agree well with 
known late-M and L dwarfs. 

A metallicity measurement for the primary star is relevant to constraining the forma- 
tion of CD-35 2722 B. Ho wever, we see that CD-35 2722 A has [Fe/H] = 0.04±0.05 dex 
( Viana Almeida et al.l 2009 ), which is close to the mean [Fe/H] of the solar neighborhood, 
—0.05 dex (j Johnson &: Appsll2009l ). As a consistency check, we make an indirect estimate 
of the metallicity of CD— 35 2722 A using the photometric calibration of I Johnson fc Apps 
( I2OO9I ). Following their method, we determined the [Fe/H] of CD— 35 2722 A from its 
height above the main-sequence isochr one in a plot of M k versus V — K. Since Ml dwarfs 
reach the main sequence by 100 Myr f lSiess et al.ll2000l ). the comparison of CD— 35 2722 A 
with nearby main-sequence stars is appropriate. We derived an [Fe/H] of 0.0±0.16 dex for 
CD— 35 2722 A, consistent with the pub l ished value. An alternate calibration of the metal- 
licity is found in ISchlaufman fc Laughliru (120 lOf ) , who use a volume-limited and kinematically 
matched sample of F and G dwarfs from the Geneva-Copehnagen Survey to infer the mean 
metallicity of M dwarfs in the Solar Neighborhood and base their final calibration purely 
on high-resolution spectroscopy of FGK primaries with M dwarf companions. Using their 
method, we derived an [Fe/H] of -0.12±0.16 dex for CD— 35 2722 A, again consistent with 
a solar metallicity. 



3.3. Spectral Analysis 

3.3.1. Spectral Typing 

To derive an empirical spectral type for CD— 35 2722 B, we compared our Gemini/NIFS 
spectra to field L dwarfs from the IRTF Spectral Library ( ICushing et al.ll2005l : iRayner et al. 
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20091 ) and young L dwarfs from recent studies (Figure [7]). When comparing to lower resolu- 
tion spectra (IRTF prism R=150, IRTF SpeX R=2000), we convolved our NIFS spectrum 
with a Gaussian of appropriate width. Since the field dwarfs are old objects, we expect some 
mismatch with CD— 35 2722 B in surface gravity dependant spectral features. 

The best matched spectra we find is from the L4.5 dwa rf 2MASS J22244381 -0158521, 
which is believed to have an unusually dusty atmosphere (see [Stephens et al.ll2009l ). We note 
that almost every absorption line feature is well-matched between this dusty field object and 
CD-35 2722 B. However the NIR colors of CD-35 2722 B are a little bluer than those of 
2MASS J22244381-0158521 {E - K=0.77 vs 0.86 and J - H =OM vs 1.05; see Figure ED 
The young (f« 5 Myr) L4 planet IRXS J160929. 1-210524 b Jlafreniere et al.ll2008l . l201oh 
exhibits a narrower //-band shape than the CD— 35 2722 B spectra, but the JHK absolute 
magnitudes and J — H and H — K colors of the two objects agree completely within errors. 

The field L3 and L5 dwarfs show clear differences in line strengths, continuum shape 
and flux level in the /C-band. The CD— 35 2722 B spectrum also shows signs of youth in 
that the shape of its //-ban d spectra is more triangular than that of field L dwarfs, as 
seen in young L-type dwarfs (ILucas et al.l l2006l : iKirkpatrick et al.l l2006l : lAUers et al.l 120071 ) 
However, the //-ba nd shape is also less pe aked than the young L dwar fs IRXS J160929.1 
210524 b (~5 Mvr: Ihafreniere et al.l[201o[ ) and G196-3 B (^ 100 Myr; iRebolo et al]ll998 



Shkolnik et ahlbood l Overall, CD— 35 2722 B does not seem to have especially low surface 
gravity, since the strength of its FeH lines and K I lines are comparable to those of field 
dwarfs. 

Since the L3 and L5 spectra already differ significantly from that of CD— 35 2722 B, we 
assign CD— 35 2722 B a spectral type of L4±l. We also derive consistent spectral types of 



L2.5, L3.5, L3.5 a nd L4 using spectra. 
Allers et al.ll2007[ ). H2OA (0.56±0.03; 
2003h . and CH4 2.2/im (0.946±0.005; 



indices from the li terature: H2O 1.53)u m (1.27=b0.05; 



Geballe et al. 



McLean et al 1120031) . H2OC f0.65±0.03: McLean et al. 



20021 ) . respectively. 



We also tried to find the best matched spectra to CD— 35 272 2 B from the > 400 M, 
L, an d T dwarfs in the online archive of SpeX Prism data (e.g. iBurgasser fc McElwain 



20061 ). by scaling relative fiux levels to minimize y^. To match the low-resolution SpeX 
spectra, our NIFS spectra are first Gaussian smoothed. At this lower resolution, most 
features of the spectra are no longer visible, and so fits rely mainly on the shape of the 
continuum. Nevertheless, we observe a clear minimum in the value for spectral types 
between L3 and L6. We verified the fits by eye to check that the minima were sensible. The 
closest match es to our smoothed sp ectra, in order of increasing xi are 2 MASS J20025073 - 
0521524 (L6, IBurgasser et"aD l2008h . SDSS J232804.58-103845.7 (L3.5, Ichiu et al.l[2006[ ). 
2MASSW J0205034+125142 (L5, Reid et al. 2006) and 2MASS J21512543 -2441000 (L3, 
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Burgasser et al.ll2008l ). 



With a spectral type estimate for CD— 35 2722 B in hand, we compare its NIR luminosi- 
ties with those of field MLT dwarfs, and other young MLT dwarfs with measured parallaxes 
(see Figure [9]). CD— 35 2722 B is over-luminous by ~ 1 mag compared to L4 field dwarfs in 
the JHK-bands. As noted earlier, its NIR luminosities are indistinguishable from those of 
the 5 Myr old pl anet, IRXS J160929 .1-210524 b. The young L4 dwarfs 130948 B and C 
(age ~ 790 Myr, iDupuy et al.ll2009al ) are also over-luminous by ~ 0.5 mags. Young late-L 
and early -T companions are clear l y und er-lur ninous compa r ed to eld objects, as previously 
noted by iMetchev fc Hillenbrandl (120061 ) and iBowler et al.l (120101 ) , while the young late-M 
to mid-L compan ions are overluminous. The young LI dwarf, AB Pic b (age ~ 30 Myr, 
Torres et al.ll2008l ) seems under- luminous and thus unlike other young M7-L5 dwarfs. 



3.3.2. Model Atmosphere Fitting 



To constrain the atmospheric properties of CD— 35 2722 B, we compare our NIFS spec- 
trum to the Ames-Dusty synthetic spectra of lAUard et al.l (120011 ). For each model spectrum 
in the Dusty grid, the flux of the model is adjusted to minimize the value of the fit to the 
observed spectrum. The overall minimum value across all template spectra in the grid 
gives the best-fit model (Figure [HD. We find the smallest values for Te// = 1700-1900 K 
with log(g) = 4.5-6.0. Since, Ames-Dusty models are not able to reproduce very well the 
spectra of cooler objects near the L-T transition, the lower temper ature bound should not b e 
taken as a very stron g constraint. Nevertheless, using models from lSaumon &: Marleyl (120081 ). 
Cushing et al.l ( l2008l ) are able to fit the spectra of L4 and L5 dwarfs with their Teg =1700 K, 
log(g)=4.5 models and L6 dwarfs with their Tegf=1400 K, log(g)=4.5 atmospheric models, 
achieving good-fits past spectral type L8. This suggests that our lower bound of Teff=1700 K 
is reasonable. Our temperature estimate is also consistent with the temperature we derive 
from the Lyon/Dusty evolutionary models in § 13.41 

Since values depend more on the quality of matches between the shapes of the 
continuum rather than matches between specific spectral features, we also judge the best-fit 
models by eye. Examining the fits, we note that while no single model is a perfect match to 
the CD— 35 2722 B spectra, the 1700K, log(g) = 4.5 model reproduces the overall continuum 
shape from J to K rather well (Figure |H]). At the same time, the shape and the depth of 
the two sets of K I doublets at 1.17 and 1.25 fim and the FeH lines in the J-band are clearly 
better reproduced in the 1700 K model than in other models. Moreover, given the young 
age of the system, we expect it to have a low surface gravity. The lower gravity of the 
1700 K, log(g) = 4.5 model is also consistent with the value estimated from the Lyon/Dusty 
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evolutionary model using our bolometric luminosity and age estimates (§ 13. 4p . Since the 
quality of the J-band fits are clearly much worse for the 1700 K, log(g)=4.0 and log(g)=5.0 
models, we estimate a log(g)=4.5±0.5 for CD— 35 2722 B. 



We also computed xt from the fits to L dwarf models with clouds presented in lBurrows et al. 



(120061 ) . The shapes of the continua and the relative fluxes in the JHK-hands are completely 
consistent between these models and the CD— 35 2722 B spectra. The strengths of the K I 
doublets in the J-band are also consistent. The xt results show a clear minima in the T^s- 
log(g) space, and we estimate Tg// =1700±100 K, and log (g) =4.5-5.0 from the Burrows 
models. The iJ-band peaks in the models are, however, clearly rounder with more flux at 
the blue end compared to the CD— 35 2722 B spectra due to the absence of FeH features in 
the models. The models also show stronger water absorption at the red end of the H-hand. 
Moreover, the models are deficient in Fe and FeH features in the J-Band, and th e ir CO 
features in the K-hand are too strong. The cloud-free models from iBurrows et al.l (120061 ) 
have very different continua shape from the CD— 35 2722 B spectra, and a reasonable match 
was not found. 



3.3.3. Radial Velocities 

Due to the dither pattern in the first set of J-band observations, there was enough light 
in the NIFS data to extract a spectrum of CD— 35 2722 A. We used this spectrum to test 
whether CD— 35 2722 A and B share the same radial velocity. The spectral types of the 
two objects are sufficiently different that cross-correlation techniques are not easily applied. 
Instead, we fit Gaussians to four strong K I lines in the J band, taking the average of the 
centroids as our measure of radial velocity. We find that the velocities of A and B agree 
to within the errors (±10 km/s) further confirming their physical association. The system 



veloc ity (25±10 km/s) also agrees with the published value (31.4±0.4 km/s, iTorres et al. 



2006! ) and the predicted value in § 3.1 within our errors. 



3.4. Mass Estimates 



We first determined the mass of CD— 35 2722 B using the Lyon/Dusty evolutionary 
models ( IChabrier et al.ll2000l ). based on the observed absolute magnitudes and estimated 
age. To compute the mass and its associated uncertainty, we used a Monte Carlo technique, 
drawing the input values for Mj and age from random distributions. Values for Mj were cho- 
sen from a normal distribution corresponding to the photometry in Table [U Age values were 
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chosen from a uni form distribution between 50 and 150 Myrs, as estimated for the AB Do- 
radus association ( IZuckerman et al.ll2004t iLuhman et aLll2005l : IClose et al.ll2007l ). Using the 
Lyon/Dusty evolutionary models in the MKO photometric system (I. Baraffe, private com- 
munication), we computed the mass of CD— 35 2722 B from 10,000 random trials, resulting 
in an adopted mass of 31±7 Mjup from the median and standard deviation of the trials. 
Using the same approach with our Mh and M^^measurements, we estimate masses of 30±7 
and 34±9 Mj„j,,respectively. Thus, CD— 35 2722 B is probably a brown dwarf, if we follow 
the notional dividing line between planets and brown dwarfs set by the deuterium-burning 
limit of 13.6 Mjup for solar metallicity. 

To assess the uncertainty contributed by each input parameter (age, distance, and J 
mag) in our mass calculation, we also computed the mass of CD— 35 2722 B by assuming two 
of the inputs were error-free. The mass uncertainty arising from the assumed age dominates 
over that from the distance, which in turn dominates over that from the J-band photometry 
(Table H]). Of course, this does not account for the systematic errors in the evolut i onary 
models thernselves, which are just now being rigorously tested (e.g. IStassun et al.l 120071 : 
Dupuv et allboOQbh . 



A somewhat more robust way to derive the mass of CD— 35 2722 B is to use its bolomet- 
ric luminosity and age as inputs to the Lyon/Dusty evolutionary models, since bolometric 
luminosities are less subsceptible to uncertainties in the model atmospheres than single band- 
pass fluxes (jChabrier et al.ll2000l ). Using the J-band bolometric corrections from Liu et al. 
(2010) for spectral types of L3-L5 (§ 13. 3p and the input distribution of Mj described above, 
we calculated the corresponding distribution of bolometric luminosities for CD— 35 2722 B. 
Then, using the luminosity and age distributions with the evolutionary models, we obtained 
a mass of 31±8Mj„p, again using the median and standard deviation of the trials (Figure [T0|) . 
The Lhoi-dehved masses from the H and K-bands were also 31±8Mjup, which is expected 
since most of the mass uncertainty arises from the uncertainty in age. We adopt this L^oi- 
derived mass as our best estimate in Table [T|, though the mass estimates from individual 
filters and from Lboi are completely consistent. Note that the same set of input data and 
evolutionary models give an effective temperature of 1980±100 K and a surface gravity with 
log{g) = 4.54 ± 0.13 dex, in reasonable agreement with the results from model atmosphere 
fits to the near-IR spectra (§ 13. 3p . 

Using the same methods as for the L yon/Dusty models, we coni pare the age and lumi- 
nosity of CD— 35 2722 B to the models of iBurrows et al.l (jl997l . l200ll ) to estimate the mass, 
Teff, and log{g) to be 36 ±6Mj„p, 2150±90 K and 4.83±0.13 dex, respectively. 

An alternative means using evolutionary models to estimate the mass comes from using 
the age and effective temperature, with the latter estimated from model atmosphere fitting 
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of the NIR spectra (§ 13. 3p . We adopt a similar Monte Carlo approach, assuming a uniform 
distribution in T^s from 1700-1900 K, and derive a mass of 23 ± 5 Mj^p- This is somewhat 
lower than the mass derived from {L}joh age}. This discrepancy occurs because the model 
atmosphere-derived T^^ is ~ 180 K cooler than the evolutionary model-derived T^^ above, 
leading to a lower massEI 

A third means to estimate the mass from the evolutionary models comes from the H-R 
diagram, namely using the bolometric luminosity and the estimated temperature. This is 
the least useful method, as small uncertainties in T^^ translate into large mass uncertainties, 
given the tightness of the model tracks on the H-R diagram. Assuming a normal distribution 
in luminosity and a uniform distribution in Tefj, we derive 11 ± 7 Mjup- Figure [TO] illustrates 
the observed properties of CD— 35 2722 B relative to the model predictions for different pairs 
of observable properties. 

The color-magnitude diagrams in F igure [TT] compare t he 50 and 150 Myr isochrones 



from the Lyon/Dusty and Lyon/COND (iBaraffe et al.l 120031 ) evolutionary models with the 



color and absolute magnitudes of CD— 35 2722 B. In both plots, the Mj vs J — H and the 
Mk vs J — K diagram, the companion's position is within Icr of that of 22, 31 and 39 Mj^p 
Dusty model-simulated brown dwarfs at ages 50, 100 and 150 Myrs, respectively. 



4. Conclusions 

As part of the Gemini NICI Planet-Finding Campaign, we have discovered a substellar 
companion to the young MIV star CD— 35 2722 A, a member of the AB Doradus association. 
High quality astrometry over one year confirms CD— 35 2722 B as a bound companion, with 
a projected separation of 3.17" (67±4 AU) in January 2009. We estimate a companion mass 
of 31 ±8 Mjup by using its bolometric magnitude and age and the Lyon/Dusty evolutionary 
models. JHK photometry show that the system has L dwarf-like colors. 

The 1.1-2.4 /.im spectra of CD— 35 2722 B reveal a low temperature and low surface 
gravity, as expected for the young age and low mass of the companion. The if-band contin- 
uum shows the triangular shape characteristic of young ultracool dwarfs, though intermediate 
in shape between the youngest known objects and old (field) objects. Through direct com- 
parison with known objects, we derive a spectral type of L4±l, making CD— 35 2722 B one 



somewhat similar effect has been identified in the young (ss 790 Myr) L4+L4 dwarf binary 
HD 130948BC, where the T^s derived from evolutionary models using the age and luminosity are « 150- 
300 K hotter than expected from atmospheric model fitting for objec ts of similar spectral type. (A direct 



model fit to the HD 130948BC spectrum is not available.) See § 5.8 of lDupuv et al 



2009a. 
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of the coolest young substellar companions identified to date. We fit the spectra using the 
Ames-Dusty model atmospheres and derive Tg// =1700-1900 K and log{g) =4.5±0.5 dex. 
These estimates are consistent with the Lyon/Dusty evolutionary model-derived tempera- 
ture of 1980±100 K and log{g) =4.54±0.14 dex. As a member of a moving group with 
a good age estimate, CD— 35 2722 B serves as a valuable benchmark for calibrating age- 
dependent spectral effects in L dwarfs. Comparison with other young substellar companions 
shows that when using near-IR spectral types, the late-M to mid-L companions appear to 
be over-luminous compared to the field objects, while the late-L and early-T companions 
appear to be under-luminous. 

The mass and separation of CD— 35 2722 B are comparable to other known young 
star-|-substellar companion systems. TWA 5b (20 Mj^p, 98 AU; Lowrance 1999), GQ Lup B 
(17 Mj^p, 100 AU; Neuhauser et al. 2005) and Gl 229 B (35 Mj^p, 45 AU; Nakajima et al. 
1995) all have masses and separations within 50% of that of CD— 35 2722 B, and all have 
primaries of spectral type K7 or later. This similarity may be a hint that they share a com- 
mon formation mechanism. In-situ formation by core accretion is not viable for such massive 
companions at such large separations, especially around an M star. Alternate theories do 
not provide easy pathways to the formation of substellar companions like CD— 35 2722 B; 
for instance, formation by fragmentation of a prestellar core (Bate 2003) is difficult for the 
mass ratio and separation of CD— 35 2722 AB. Formation via gravitational instability (Boss 
1997, Rafikov 2005, Boss 2006, Stamatellos & Whitworth 2008) require large massive disks 
(>0.2-0.35 Mq) and thus also do not provide an obvious solution, given the low mass of 
the primary star. However, Jeans mass fragmentation of the interstellar molecular cloud 
may be able to ex plain s tellar- substellar systems, e.g., as discussed by Zuckerman and Song 



(2009). Moreover, iBosd ( 1200 ll ) found that dense molecular cloud cores with magnetic fields 
can collapse and fragment into multiple protostar systems with initial masses as low a Jupiter 
mass, i.e., even smaller than the companion to CD— 35 2722 A. This suggests that the CD- 
35 2722 system might well have formed by the classic binary star formation mechanism of 
fragmentation of a collapsing dense molecular cloud core. 
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Fig. 1. — A cutout from the reduced NICI image of CD-35 2722, representing a 45- minute 
ADI observing sequence (North is up, East is left). This image is the sum of the reduced 
images from the blue and red channels (CH4 4% short and long filters). The stellar position 
is indicated by a 'star' symbol. CD— 35 2722 B was detected 3.137" away from the primary 
at a PA of 243.10° on UT Jan 10, 2010. The background object, shown in the square, had 
5.316" separation and PA=276.32°. The pixel intensities are shown in linear stretch from 
— 4cr to 4cr, where a is the pixel-to-pixel noise in background. 
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Fig. 2.— NICI images of CD-35 2722, obtained on UT March 9, 2009 in the J, H and 
Ks-hsmd filters. North is up, East is left. The primary star, attenuated by the partially 
transmissive coronagraphic mask, is in the upper left area of the panels. CD— 35 2722 B is 
in the lower right area, roughly 3" away from the primary. 



Fig. 3.— A NIFS A-B image of CD-35 2722 B with a 3.2" x 3.2" FOV. North is up, East 
is left. The primary, not shown in the image, hes above and to the left of this region. 




Fig. 4.— Normalized JHK-hand NIFS spectra of CD-35 2722 Bfred), plotte d along with 
the dusty L4.5 dwarf 2MASS J2224438 1-0158521 (blue) knapp et al.l l2004h. Molecular 
features and line identifications are from ICushing et al.l (120051 ). 
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Fig. 5. — Comparisons of the observed motions of the two candidates (Left: CD— 35 2722 B; 
Right: background object) to the expected motion of a background object given the first 
epoch astrometry, the proper motion and the distance of CD— 35 2722 A. The very small ob- 
served change of the PA of CD— 35 2722 B (bottom left panel) differs significantly (3cr) from 
that expected for a background object. The change in the position of the other candidate is 
consistent with a background object. 
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Fig. 6. — The infrared colors of CD— 35 2722 B (center red square) plotted along with the 
colors of field objects and low mass companions. The JHK colors of CD— 35 2722 B are 
consistent with a late-type dwarf, but slightly bluer than that of 2MASS J22244381-0158521 
(dark blue triangle), the dusty L4.5 with very similar spectra. The colors also agree com- 
pletely with those of IRXS J160929. 1—210524 b, which is plotted as a grey rectangle, almost 
covered by the symbol of CD— 35 2722 B. The colo rs of the field objec ts are derived from 
synthesized photometry of publis hed spectra (from iRavner et al.ll2009l ). Th e L dwarf col- 
o rs come from the compilation bv iLeggett et al. (based on data from iGeballe et al.l 



20021 : iGohmowski et al.ll2004l : iKnapp et al.ll2004J : IChiu et al.ll2006f ). with the plotted symbols 
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Fig. 7. Gemini NIFS spectra of CD-35 2722 B (black curves) compared to young (low gravity or dusty) and field L dwrarfs. The J, H and /C-band spectra 

are sfiown in tfie left, middle and right panels, respectively. The NIFS spectrum is smoothed to match the resolution of the comparison spectra, which have 
been scaled to fit each band separately. The CD— 35 2722 B spectrum is most similar to that of the dusty L4.5 2MASS J22244381-0158521 and the low gravity 
L3 dwarf G196-3B, and begins to differ significantly from spectra which are more than one subclass away from L4. Thus we assign a spectral type of L4±l 
to CD-35 2722 B. The top five comparison spectra (2MASS J15150083+4847416, SDSS ,1053951.99-005902.0, 2MASS J15074769-1627386, 2MASS J22244381- 
0158 521, 2MASS J00361617 +1821104) are from the IRTF spectral library (R = 2000 ; [Custiing et al. 2005). The IRXS J160929. 1-210524 b spectra (R = 850) are 
from lLafreniere et al.l ||2008| ). while the G196-3 B spectra are from lAUers et al.l l l201(t) . Molecular features and line identifications are from lCushing et al.l bOO^ . 
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Fig. 8. — Left: The xZ from fitting Ames-Dusty model atmospheres (lAllard et al.ll200ll ) to 
the NIFS spectra of CD— 35 2722 B as a function of model temperature and gravity. The 
best fit models are the Tg// = 1700-1900 K models with log(5r) = 4.5 — 6.0 dex. However, 
none of the models match the continuum shape very well. The reduced chi-squared values 
quickly deteriorate for Teff < 1600 K and Teff > 2000 K models. Right: The best fit 
Ames-Dusty model spectra for different effective temperatures. 
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Fig. 9. — The absolute J, H and K MKO magnitudes and spectral type of CD— 35 2722 B 
compared to cool field objects and other young MLT dwarfs with parallax measure- 
ments. In the NIR, CD— 35 2722 B is over-luminous compared to field objects, but has the 
same JHK luminosities as the 5 Myr old Upper Sco planet, IRXS J160929.1-210524 b. 
Earlier young dwarfs are generally also over-luminous, while later youn g dwarfs are 



under - luminous. The objects shown are the MLT dwarf compil a tion f rom iLeggett et al. 
(120101 ) (gray symbols). Upper Sco L dwarfs f rom Lodieu et all (120081) fgreen symbols) 
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Fig. 10. — Comparison of the properties of CD— 35 2722 B to those predicted by the Dusty 
model mass tracks. The mass estimate derived from the bolometric luminosity and age 
and from the Tg// (derived from model atmosphere fits to the near-IR spectra) and age are 
consistent with each other. 
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Fig. 11. — Comparison of the COND and Dusty model isochrones at 50 Myr, 150 Myr 
and 1 Gyr with the absolute magnitudes and NIR colors of CD— 35 2722 B shown by the 
red square. For comparison, we also plot the NIR magnitudes and colors of field M and 
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not predict the field brown dwarf sequence well at low temperatures. 
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Table 1. Properties of the CD-35 2722 AB system. 



Property CD-35 2722 A CD-35 2722 B 



Age (Myr)^ 


100±50 


Distance (pc)'' 


21.3±1.4 


Angular separation ("Y 


3.137±0.005 


Physical separation (AU)*^ 


67±4 




Position angle {°Y 


243.10±0.25 


A J, IRTF (mag) 


5.77±0.06 


AH, IRTF (mag) 


5.50±0.06 


AK, IRTF (mag) 


4.98±0.05 


A J, NICI (mag) 


5.66±0.12 


AH, NICI (mag) 


5.48±0.15 


AKs, NICI (mag) 


4.96±0.07 


AK, NICI (mag)'^ 


4.94±0.07 


Jmko (mag)'' 


7.86±0.03 


13.63±0.11 


Hmko (mag)'' 


7.28±0.03 


12.78±0.12 


Kmko (mag)'' 


7.03±0.05 


12.01±0.07 


Mj (mag) 


6.0±0.14 


11.99±0.18 


Mh (mag) 


5.4±0.14 


11.14±0.19 


Mk (mag) 


5.1±0.14 


10.37±0.16 


(J - H)mko (mag) 


0.58±0.03 


0.86±0.13 


{H - K)mko (mag) 


0.25±0.02 


0.77±0.2 


Spectral Type 


MlVef 


L4±l 


Mm (mag) 


8.11±0.25 


13.7±0.3 s 


Mass 


0.4±0.05 Mq^ 


31±8 Mjup' 


Te// (K) 


3680±100 j 


1700-1900 ^ 


log(fi') (dex) 




4.5±0.5 ' 


[Fe/H] 


0.04±0.05 dex 





Age estimat e for the AB D o radus association 



(Zuc 


ierman et al. 


2007 


)• 



2004; Luhman et al. 2005; Close et al. 
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^Distance from parallax measurement using CAPScam. See 
§3.1 for details. 

"Epoch January 10, 2010 UT. 

"^Tra nsformed from NICI magnitude, using the relation- 

ship in [Stephens fc Leggettl (l2004j ) and asumming a spectral type 
of L4. 

*^MKO magnitudes de rived from 2MASS magnitudes using 
color transformations in iLeggett et al.l ( 120061 ). All photometry 
and colors are given in the MKO system. Apparent and absolute 
magnitudes, colors, and mass estimates for CD— 35 2722 B are 
all derived from IRTF photometry. 



^Spectral Type from [Torres et al.l ( 20061 ). 



^Bolometric magnitude of primary derived using bol o metri c 
magnitude equations for M dwarfs in ICasagrande et al.l ( l2008l ). 
Bolometric magnitude for companion derived from Mj, based on 
bolometric corrections in Liu et al. (2010). 



^Mass estimate based on lSiess et al.l ( 12000! ) using Mh and age. 
'Mass estimate based on Lyon/Dusty model using Mfeo/ and 



age. 



JTe// for an MIV from lGautier et al.l (120071 ) 



'^Tg// from Ames-Dusty atmospheric model fitting result. 

' Log(g) from fitting Ames-Dusty model atmospheres to the 
Gemini/NIFS near-IR spectra. 



[Fe/H] measured by lViana Almeida et al.l (l2009l ) 
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Table 2. Astrometry of objects near CD— 35 2722 A. 



Object 


UT Date (Y/M/D) 


Separation (") 


Position Angle (°) 


CD-35 2722 B 


2009/01/06 


3.172±0.005 


244.13±0.25 




2010/01/10 


3.137±0.005 


243.10±0.25 


background object 


2009/01/06 


5.312±0.011 


275.48±0.25 




2010/01/10 


5.316±0.011 


276.32±0.25 



Table 3. Spectroscopic observations obtained with Gemini- North/NIFS. 



UT Date (Y/M/D) 


Band 


No. of Exp. 


Exp. Time(s) 


Comments 


2010/03/07 


J 


12 


300 




2010/03/18 


H 


9 


300 


Partial cloud 




K 


4 


120 




2010/03/19 


J 


9 


300 


Poor conditions, partial cloud 



Table 4. Comparison of the sources of uncertainty in CD— 35 2722 B's estimated mass. 



Parameter 


Measurement 


Derived Mass Error 


Age 


100±50 Myr 




Distance 


21.3±1.4 pc 


0.6 Mj^p 


Apparent J mag 


13.63±0.11 mags 


0.3 Mj^p 


All 




8 Mjup 



Note. — For each parameter, we set the uncertainty in 
all the other parameters to zero and calculated the error 
in CD— 35 2722 B's mass using the Lyon/Dusty evolutionary 
models. 



